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I. INTRODUCTION 

When a system described by an order parameter, ini- 
tially placed into a high temperature single phase region 
of its phase diagram, is brought to a point inside the co- 
existence curve by a sudden change of temperature it be- 
comes thermodynamically unstable and phase separates 
as a result of the existence of many competing ground 
states. After the quench the system can order kinetically 
through either nucleation or spinodal decomposition. In 
the latter process a microscopic long-wavelength fluctua- 
tion initially present is amplified and determines the for- 
mation and evolution of various patterns characterized 
by the presence of a universal length scale L(t), associ- 
ated with the typical domain-size and separation among 
topological defects. As the system orders, L(t) grows in 
time in a power-law fashion t 1 ^ and the time-dependent 
structure factor C(k,t) displays dynamical scaling. 

A successful approach to the study of these phenomena 
is represented by the time dependent Ginzburg-Landau 
equation. Many years ago Hohenberg and Halperin 
provided a useful classification of the various models, 
which comprises a vast class of dynamic critical phenom- 
ena, in terms of few parameters. Within their scheme two 
models have received a great deal of attention: model A, 
where a single field evolves towards equilibrium with non 
conserved order parameter dynamics (NCOP) and model 
B where the order parameter is conserved (COP) 

The NCOP dynamics is aimed to describe an ordering 
process similar to that of the Ising model with conven- 
tional Monte Carlo spin flip dynamics, while COP ac- 
counts for the approach to equilibrium of an alloy. In 
this case, one observes that the system orders by grow- 



ing droplets larger than a critical size at the expenses of 
smaller droplets, while keeping the total amount of ma- 
terial fixed. The effect of the conservation law is to slow 
down the phase separation process, because the material 
has to be transported via diffusion through the system 
before being added to a growing region. The value of the 
dynamical exponent z is z — 2 for NCOP, whereas for 
conserved dynamics is z = 3 for a scalar order parameter 
and z = 4 for vector order parameter, indicating that 
conservation laws play an important role in the dynami- 
cal process. 

While COP and NCOP dynamics has been widely 
studied, the Phase field model ||, which describes the 
coupling of a NCOP system with a diffusive COP field 
such as temperature or concentration, seems not to have 
been completely explored, in spite of the fact that it dis- 
plays a variety of interesting peculiar features. Only 
to mention the most striking of these, we recall that 
the Phase field model accounts for the regularity of the 
shapes observed during the growth of crystals into a su- 
percooled melt. According to the Phase field model a pla- 
nar solid front growing in the supercooled liquid under- 
goes the so called Mullins-Sekerka instability. This phe- 
nomenon can be understood as follows. The latent heat, 
released when the liquid freezes, is diffused into the colder 
liquid and thus promotes the freezing of more material. 
The larger is the temperature gradient the faster is the 
advancement of the front. Now, imagine to slightly per- 
turb the isothermal flat solid-liquid interface in a slowly 
varying fashion. As a result of the deformation, the tem- 
perature gradient will be larger on the bulges of such a 
boundary and so the heat flux. This fact makes a solid 
tip to grow faster than a flat portion of interface and 
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provides a mechanism by which a perturbation of finite 
wavelength is amplified, as it was discovered by Mullins 
and Sekerka g§. 

In a recent letter || we have introduced an N- 
component version of the Phase field model in order to 
study the evolution of a non conserved order parame- 
ter bilinear ly coupled to a conserved field U. In the 
limit N — ■+ oo we were able to obtain some analytical re- 
sults on the non-equilibrium relaxation behaviour of the 
model. Here we expand these results. 

The model, inspired by the model C 0, displays sev- 
eral interesting features: the evolution of the vector field 
4> is non conserved in the early regime, then after a 
crossover time it develops an instability at finite wave- 
length due to the coupling with the conserved field U. 
In the very late regime the COP behaviour becomes 
eventually dominant and <fi shows a genuine COP evo- 
lution, including multiscaling. A similar mechanism was 
reported by Somoza and Sagui in a numerical study of 
the model C where they observed that notwithstanding 
the non conserved field evolves faster than the conserved 
field. For late times the growth is driven by diffusion of 
the conserved variable and the order parameter becomes 
slaved by the diffusive field. 

The paper is organized as follows. In section || we give 
a brief physical motivation of the model. The equation 



defining the model are given in section [II, and its equi 



librium and dynamical properties are discussed in section 
IV and section [v[ respectively. Finally section VI con- 
tains a brief summary of the results and discussions. 



II. MOTIVATION OF THE MODEL 

The model we discuss in the present paper belongs to 
the family of the O(N) spherical models and it has been 
introduced Q with the aim of studying exactly the cou- 
pling of a NCOP field with a diffusive COP field. The 
O(N) generalization proves to be fruitful because, while 
retaining the salient features of the phenomena occurring 
during the diffusion limited growth it allows for some an- 
alytical results in the limit N — > oo. 

Historically, models containing couplings quadratic 
with respect to the material NCOP field <f> and linear 
with respect to the COP diffusive thermal field u were 
introduced as early as the seventies in the framework 
of the dynamical critical phenomena and named models 
C. Few years later, Langer, in order to study first order 
phase transitions accompanied by latent heat of fusion, 
put forward the so called Phase field model || , in which 
the coupling was assumed to be bilinear with respect to 
the two fields. 

The material is characterized by an order parameter <j) 
which assumes a positive value in the solid phase and a 
negative value in the liquid phase. The local tempera- 
ture of the system is treated as an additional dynamical 



field obeying a heat diffusion equation in the presence of 
sources represented by the amount of material changing 
phase. The solidification takes place adiabatically so that 
no heat can flow to the outside. One defines the dimen- 
sionless temperature field as u(x, t) = c p (T(x, t)—T m )/L, 
where L is the latent heat of fusion per mole, c p the spe- 
cific molar heat at constant pressure and T m is the bulk 
melting temperature. The spatial average of u at the ini- 
tial time is the so called undercooling parameter A and 
is a negative quantity. 

The thermal field diffuses according to the modified 
Fourier equation: 



du(x, t) 
dt 



DV 2 u{x,t) 



d<p(x,t) 
dt 



(1) 



where D is the thermal diffusivity and the last term on 
the r.h.s. is the amount of material which crystallizes 
per unit time and thus proportional to the heat released 
during the first order transition. Both c p and D are as- 
sumed to be equal in the two phases. The evolution of 
the order parameter is determined by the non linear time- 
dependent equation of the Ginzburg-Landau type [p|-[i"o[ : 



dcj)(x, t) 
dt~ 



-r 



F[<]>,u] 



x,t) 
r[-V 2 + r0 



(2) 



To describe a two phase system the form of F is chosen to 
be a double well for r < 0. The coupling to the thermal 
field u can create an unbalance, in such a way that for 
negative values of au the liquid phase {4> ~ — y—r/g) 
is metastable with respect to the solid {<fi ~ + y/—r/g). 
In the absence of coupling to the temperature field, i.e. 
a = 0, eq. (Q) represents the familiar Cahn-Allen equa- 
tion, also called Model A. 

The process contains two stages: during the first stage 
the solid grows at the expenses of the liquid, while in the 
second stage the total amount of solid is nearly constant 
and the growth is limited by diffusion of the thermal field. 

Interestingly, the two dynamical equations (|l|) and (||) 
can be derived from a Lyapounov functional T . which 
plays the role of the time dependent Ginzburg-Landau 
potential in the present problem. If one performs the 
transformation U = u — <f> and eliminates u in favor of 
the new field U one can write eqs. dl|)-(0) as 



dcj>(x, t) 

dU(x,t) 
dt 



-r 



(x,t) 



u 



D „ ST 
— V 2 

a SU(x,t) 



(3) 
(4) 



with the Lyapounov functional 
T\4>, U] -- 



d d x 



9 ,4 



Q 



(5) 
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The functional T has two equal minima when the tem- 
perature field vanishes, i.e., for U = —(f) and generates a 
complex dynamical behavior which has been the object 
of some studies. However, its global properties are not so 
well known. This fact, lead us to formulate an O(N) in- 
variant vectorial generalization of the above model. This 
kind of models, in fact, lend themselves to nearly an- 
alytical solutions thus providing useful insights on the 
properties of the scalar order parameter solutions. 



III. THE 0(N) MODEL 

We shall consider a system described by two coupled 
TV-component vector fields = ((/)i(x, t), <j)^{x,t)) 
and U = (Ui(x, t), Un(x, t)), whose Hamiltonian can 
be represented by fll]J|: 



H[<j>{x), U(x,t)] 



;<t> 2 



d d x 



2\2 



AN 



5(V0) 2 



2 



(6) 



where r and g, with g > and w > 0, are the standard 
quadratic and quartic couplings of the Ginzburg-Landau 
model and the last term represent a bilinear coupling be- 
tween the field <f> and U. The first three terms in eq. (j^) 
constitute the familiar Ginzburg-Landau- Wilson Hamil- 
tonian describing an O(N) <fi 4 model, whereas the last 
two terms represent the interaction between the order 
parameter field and an external fluctuating field U. 

We assume that the field <fi evolves according to NCOP 
dynamics: 



d(f>a{x,t) 

dt 



b a {x,t) 



H[tf»,U\+ri a (x,t) (7) 



whereas the field U is conserved and relaxes according 
to: 



dU a (x, t) 
dt 



rvv 2 



SU a {x, t) 



H[<f>,U]+£ a (x,t) (8) 



The noises appearing on the right hand sides of eqs. 
(B) have zero average and two-point correlations: 

(rj a (x,t)i ll 3(x',t)) = 

2T f T 4> 5 at p5(x - x')5{t-t') 

«)£/»(*',*)) = 



(9) 



(11) 



where Tt is the temperature of the final equilibrium state 
whereas Tu and are the kinetic coefficients. 

Introducing the Fourier components of the fields one 
can write the evolution equation as, 



d<j> a (k,i) 
dt 

dU a {k,t) 
dt 



= F£(k)+ Va (k,t) 



(12) 
(13) 



where F^jj are the Fourier transforms of the first term 
on the r.h.s of eqs.(fj]) and (||). 

In the limit N — > oo the cubic term entering into F^ 
can be decoupled and we have: 

F${k)=M^{ht)<l>a{k,t) + Mt U {k,t)U a {k,t) (14) 
Fg(k) = M U4> (k,t)cj) a (k,t) + Muu(k,t)U a (k,t) (15) 

where the matrix elements are given by 



Mw(k,t)= -r [fc 2 + 

M^u(k, t) = -T^fj. 

M u<t> (k,t)= -Tunk 2 , 

Muu(k,t) = -T uW k 2 . 



(16) 



The quantity S(t) is the integrated ^-structure function 
1 w 

S(t) 



Jj^2(<l>ot(x,t) (f> a (x,t)) 



d d k 
(2n) d 



b a (k,t) (f> a (-k,t)) 



(17) 



and the integral contains a phenomenological momentum 
cutoff A. The average is over the external noises r\ and 
£ and initial conditions. 

To study the behaviour at finite temperature 2/ it is 
useful to introduce the equations of motion for the three 
equal-time real space connected correlation functions 
C H (r,t) = (<t> a (R + r,t)4> a (R,t)) C^ v {r,t) = {4> a {R + 
r,t)U a (R,t)) and C uv (r,t) = (U a (R + r,t)U a (R,t)}, 
whose Fourier transforms are the structure functions. 
These correlations are independent of the index a due 
to the internal symmetry. In the N — > oo limit the struc- 
ture functions evolve according to the following set of 
equations: 

1 d 

2 ^CW(M) = M^(M)CW(M) 

+ M <t>u (k,t)C <l>u (k,t)+T ! i > T f (18) 



— C#KM) = AMM)CW(M) 



+ (Muu(k,t) + M^(M))CW(M) 
+ Afytr(M)CW(M) (19) 



-2T f T u 5 a , V 2 8{x-x')8{t-t') (10) 1 d 



Cuu{k,t) = M u4 ,(k,t)Cw(k,t) 

+ Muuik^Cuuik^+TuTjk 2 . (20) 

In what follows we shall be more interested into the 
behaviour of the field 4>, since it is the relevant order 
parameter of the system. 
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IV. EQUILIBRIUM PROPERTIES 

In this section we investigate the equilibrium proper- 
ties of the model @ and (||). It can be shown that the 
random process characterized by the Langevin equations 
(p|)-([T3|) obeys detailed balance since the following "po- 
tential conditions" , analogous of the Onsager relations, 
are fulfilled fTJ: 



5(j) a (k) 
S 

S(j) a (k) 
6 



su a{ k) F ^) 



S(pp(k 

r su p {k') 

k' 2 6 



(21) 
(22) 
(23) 



If detailed balance holds, the stationary probability den- 
sity reads 



P st [cf>,U} =Nexv\-—H[<t>,U] 



1 



(24) 



where J\f is a normalization constant. 

The equilibrium probability density is quadratic in the 
field U, therefore as far as the static properties of 4> are 
involved, the field U can be traced out. One is then left 
with an effective Hamiltonian for the field <p: 



H cB [4>] 



d d x 



(V0) 2 



-4> 2 



(0 2 ) 2 



AN 



(25) 



where r c g — r — \j? /w is the "renormalized mass". The 
importance of the field U can be fully appreciated only 
in the non-equilibrium dynamics of the system, as will be 
discussed in the next section. 

Before considering the dynamics, we briefly discuss the 
static properties of this model. Unlike the case where 
the U is quenched []l3| the system displays for space di- 
mensions d > 2 an order-disorder transition when Tf is 
lower than the critical temperature T c . In order to locate 
the critical surface Tf — T c {r, g, (J,,w) one considers the 
long-range behaviour of the structure functions C^(k), 
which can be computed from ( p5| ) . The fourth order term 
makes the calculation difficult for finite N. However for 
N — > oo we can use the Hartree approximation, exact in 
this limit, and we readily obtain for T > T c : 



phase. The critical temperature is given by the usual 
form of the 4 theory T c = (fi 2 /w - r){d-2)/(gK d ^ 2 K d ) 
with 1/K d = 2n d /' 2 T(d/2) where T(x) is the Gamma 
function. 

For temperatures Tf below T c there exists a non van- 
ishing order parameter M = (4>i), which can be assumed 
to be directed along the a — 1 direction without loss 
of generality. The (N — 1) components of the correla- 
tion function orthogonal to the order parameter direction 
diverge at small k, reflecting the existence of Nambu- 
Goldstone modes, i.e, excitations of vanishing energy cost 
in the long wavelength limit. The real space two point 
correlation function takes the form 



(0 Q (r) cj> a {r)) 



M 



fJ 



M 2 6 n 



w 



Soo(*) 



Tr 



(28) 



(29) 



where defined in equation (27) comes from the trans- 
verse components only. 

The other equilibrium correlation functions can also 
be obtained from the stationary equilibrium distribution 
and read for T > T c : 



Cuu(k) 



w — [i 2 {k? + r + gSao) 1 



cv(fc) = - - 



w' k 2 + r + gSca — fi 2 /w 



(30) 
(31) 



Note that both C^,£/(fc) and Cuu(k) are singular for 
k 0. 

We conclude by noting that to obtain the static struc- 
ture functions from the dynamical equations one has to 
supplement the requirement that the right hand sides of 
cqs.(|l8|)-(20) vanish with the following stronger condi- 
tion: 



lim 

t — >oo 



+M uu (k,t)C 4>u {k,t) =0, (32) 



C^(k) 



r + gSoc 
d d k 



fc|<A 



(2tt) 



fj? /w 



(26) 
(27) 



For Tf < T c the structure function diverges at small k 
because the full mass term r + gSoc — [i 2 /w vanishes, sig- 
naling the appearance of the ordered phase. In fact, the 
model for r < fi 2 /w and g > displays a high tempera- 
ture paramagnetic phase and a low temperature ordered 



lim M (j>lj> (k,t)C ( j, u (k,t) 



+M, pu (k,t)C uu {k 1 t) 



(33) 



to ensure that the equilibrium properties of the model 
are independent on the kinetic coefficients T^, and Tjj. 
The conditions (^2j)-(Q) can also be deduced from the 
the equilibrium properties of the model (see Appendix) . 
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V. DYNAMICAL PROPERTIES 

Since the behavior at Tf — is representative of the 
entire dynamics in the ordered phase when Tf < T c , we 
shall neglect the noise terms in the following analysis || . 

For general initial conditions the two fields are not in 
equilibrium, and we may expect that the relaxation of <fi 
is only slightly modified by the external field U . Since 
the dynamics of U is sufficiently slow compared with that 
of <p, the presence of U does not modify qualitatively the 
NCOP behaviour of <p. In particular, the size of the 
domains of </> should grow with a characteristic length 
L(t) ~ t 1 / 2 , while the maximum of the structure factor 
is located at k — and should increase in time with the 
power 

This kind of behavior persists until the domain size 
reaches the typical length associated with the field U, 
and given by the maximum of the structure function of 
U. At this stage the dynamics of <p slows down because 
the coupling with the conserved field U introduces an 
additional constraint on the dynamics of (p. For longer 
times the two fields equilibrate and the COP behaviour 
eventually becomes dominant. 

A simple analysis of the equation of motion for N —> oo 
gives the scaling of the crossover time with the coupling 
constant /i. Indeed it is simple to see that making the 
rescalings 



U/n -> u- 



k/ n —> k: 
r/n -» r; 
A/fi — > A; 



(34) 



the parameter /1 disappears from equations of motion for 
d> = d) a and U = U a : 



di 



<f>{k,t) 



T^k 2 +r + gS(t)](f>{k,t) 

-T^U(k,t) (35) 
-T uf ik 2 cj){k,t) -T uW k 2 U(k,t). (36) 



As a consequence the crossover time scales as l//i 2 . From 
this analysis it follows that if the dynamics of U is suffi- 
ciently slow then for 1 Ct< 1 /^i 2 the field cf> exhibits a 
NCOP behaviour while for t ^> l/fi 2 a COP behaviour. 
If the dynamics of U becomes too fast the first NCOP 
behaviour shrinks and becomes hardly observable. 

This scenario can be confirmed by solving the equation 
of motion ([35]) and ([36]) in a quasilinear approximation. 
To this end we assume that 



R(t) = r + gS(t) 



(37) 



is slowly varying in time, so that it can be considered con- 
stant over successive intervals of time. In other words, we 
make a piecewise linearization of the equation of motion 



along the trajectory. In spite of that, the approxima- 
tion is sufficient to identify the different regimes of the 
relaxation process. 

If we neglect the time dependence of R(t) and assume 
it to be nearly constant eqs. (|3f]) and ([36]) become a 
linear system whose solution has the form: 



4>(k, t) = c+(fc)e"+W* 
U(k,t) =4(fe)e w +W* 



•cI(/c)e t "-W* 
c7r \k) e"-( fe )* 



(38) 



where u>+(k) and w_(fc) are the eigenvalues of the M 
matrix, 

1 r 



T^k 2 + K)-T u wk I 



uj±(k) 



± yF4k 2 + R)+ r v wk 2 ] 2 + AT 4> T uf i 2 k 2 

(39) 

For time f > 1 the dynamical behavior of the solution is 
determined by the larger eigenvalue u + (k). For large k 2 , 
the eigenvalue uj + (k) decreases proportionally to —k 2 and 
hence large momenta are exponentially damped. More- 
over we see that ui+(k) is a function of k 2 which either has 
an extremum at k = or a single maximum for k =/= 0, 
as one can verify by inspecting the small-fc behaviour of 
cu+(k). The behavior of uj+(k) is shown in Fig. [l] for 
t < Tf and t > Tf. The crossover time r/ is defined 
as the time when the fastest growing mode moves from 
k = to k ^ 0. Other definitions of Tf are possible, e.g., 
the time when the peak at k ^ becomes higher then 
the k = one. However, all definitions lead to similar 
results. 

Below the critical temperature T c and in the early stage 
of the ordering process the value of gS(t) is small com- 
pared with r, i.e. R < 0, and the larger eigenvalue is well 
approximated by: 



.(*) =r^\R\ - r -r 



tr 

A brief calculation reveals that 
(i 



R\ 



k 2 + 0(k 4 ) (40) 



4>(k,o) 
t a r 



R 



U(k,0) + O(k 2 



I' 



#M) + ^tf(k,o) 



k 2 + 0(k 4 ). 



(41) 
(42) 



Therefore, assuming </>(fc,0) + (/i/R)U(k,0) — 0(1) for 
k — » 0, the coefficient cf(k) is finite for k — * while Cy(k) 
vanishes, indicating that the amplitudes of the longest 
wavelength components of U are decreased due to the 
conservation law. 

As a consequence, below T c the structure factor 
C0^(fc,t) develops a peak centered at k = 0, growing in 
time as a power of t. The structure functions C ( j,u(k,t) 
and Cuu(k,t) also develop a peak, centered at a finite 
value of k, say kf, as a result of the competing ef- 
fect between the k dependence of the exponential fac- 
tor exp[w + (fc) t] and the amplitude c^(fc), see eq (f4^). 



5 



This mechanism selects a set of exponentially growing 
[/-modes with wavevectors in a certain range centered 
around k / , whose dependence on the coupling /i is shown 
in Fig. |^. Such modes represent an inhomogeneity of the 
U field which in turn affects the spatial properties of the 
4> subsystem. One witnesses a strong feedback process 
between the two fields and the outcome is the slaving of 
the NCOP dynamics of the field <j> to the COP dynamics 
of the [/-field. 

The power law growth of C$$(k — 0, t) can be ex- 
tracted from the quasilinear approximation by using (40). 
In the early regime R starts from a negative value and 
grows towards zero due to the growing of C^(fc, t) for 
small k. This in turn implies that S(t) tends to a fi- 
nite value for increasing time. By imposing this condi- 
tion, and making use of ( [40| ) it follows that C^(k — 
0,t) ~ t d / 2 , as in the pure NCOP, i.e. the longest wave- 
length fluctuations grow faster. We note that while the 
quasilinear approximation leads to the correct scaling of 
C^(k = 0, t) ~ t d / 2 and of the domain size L(t) ~ t 1 / 2 , 
it gives for R(t) ~ log(f)/f which reveals that the approx- 
imation is slightly crude. 

These results are valid for R//i 2 not too large, i.e., far 
from the crossover region where R changes sign. Unlike 
the pure NCOP, where R(t) goes to zero, after a char- 
acteristic time Tf — 0(1/ /i 2 ) the value of \R\ becomes 
0(fi 2 ) and the NCOP behaviour ends. By inspection of 
cq.(^0|), we see that if \R\ffJ 2 < Tu /T^ the maximum of 
u>+(k) moves away from k = and the system looses its 
NCOP behaviour. 

This regime corresponds in our model to the instabil- 
ity which is observed in systems where a non-conserved 
order parameter is coupled to a conserved field. We must 
stress that in order to observe the NCOP behaviour the 
dynamics of U needs to be sufficiently slow with respect 
to the characteristic time Tf ~ 1/fJ 2 , whose dependence 
on \x is shown in Fig. |L 

For times t — 0(1/ /i 2 ) the quantity R changes sign be- 
coming positive, and finally for t — » oo tends to a finite 
value fx 2 /w while the maximum of u>+(k) moves again to- 
wards vanishing wavevectors. The dynamics is therefore 
dominated in the regime t ^> l/^ 2 by long wavelength 
fluctuations. We can then expand to+(k) in powers of k 
obtaining 



gs 



w I k — c^k 



(43) 



where C4 is a positive coefficient having a finite limit for 



By imposing that S(t) has a finite non zero, 



limit for t — ► 00 and making use of eq.(|43|) one obtains 
that in this regime 



C M (k,t) = [L(t) 2 k m (t) 



2 _ Wfc/fcm(t)) 



(44) 



L(t) ~ t 1 / 4 , 



k m (t) 



d \ogt 
4 ~T~ 



1/4 



(45) 



a behaviour typical of COP dynamics The multi- 

scaling function ip(x) is given by 



ip(x) = 1 ~ (x 2 - if 



(46) 



The COP behaviour is also observed if one considers the 
structure functions C(/,u(k, t) and Cuu(k, t). Such a mul- 
tiscaling behavior follows from the competition of two 
marginally distinct lengths, namely the domain size L(t) 
and fe~ . 

Finally we note that the quasilinear approximation in 
this regime leads to 



R(t) - fi 2 /w 



log* 
t 



1/2 



(47) 



where 



From equation (g6|) we see that R — /i 2 /w plays the role 
of the mass term r + g S(t) in pure COP dynamics [03, 
therefore in spite of the fact that the quasilinear approx- 
imation is quite crude, it gives nevertheless the correct 
scaling behaviour of COP dynamics. 

The above theoretical predictions were checked by in- 
tegrating numerically the system of equations (|l§|)-(|20|) 
by the Euler method. The k integrals were evaluated 
by a Simpson rule discretizing the wavevectors in the 
interval [0, A]. Figure^ displays the structure function 
C<p<p(k,t) for various values of the time t. One clearly 
sees that in the early regime the fastest growing modes 
are centered about k = 0, because long-wavelengths fluc- 
tuations of the field (f> increase more rapidly than shorter 
ones, whereas for t > Tf a finite wavevector peak ap- 
pears. Moreover, the growth, in this late regime, has a 
conserved character because its value at k — remains 
constant. The evolution of C^u(k,t) and Cuu(k,t) is 
shown in Figs. || and |[ Initially the fields U and (f> 
evolve as if they were nearly independent and the Cjju 
correlations display the usual finite wavevector peak of 
the conserved dynamics, whereas the (f> field evolves ac- 
cording to a faster non conserved dynamics. The long 
wavelength fluctuations of the U field are hindered by the 
conservation law and the presence of the C^u term has 
only a small effect on the C^. However, as the domain 
size L(t) reaches a critical value and becomes compara- 
ble with A/, the typical length of the oscillations of the 
diffusive field, the two fields strongly interact. Within 
this late regime the dynamics becomes controlled by the 
conservation law induced by the U field. In Fig. ^ we 
show the behavior of the height of the peak of C<^ versus 
time, where one clearly sees the crossover from the early 
time behavior t d / 2 to the late stage slope t d ^ 4 . In the 
crossover regime due to the presence of a double peak 
the maximum height decreases until the peak at k = 
disappears. 
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Finally we report the numerical result concerning the 
multiscaling, observed in the late regime. In Fig. || we 
display the shape function F(x) = k^ l (t)C^(xk m , t) as a 
function of x . In Fig. ^ we show the multiscaling function 
tp(x) obtained from the best fit of C^(k,t) as a func- 
tion of L(t) 2 k m (t) 2 ~ d for fixed values of x = k/k m (t). 
Similar curves can be extracted from the other structure 
functions. We note that while the data follow quite well 
the theoretical result ( f46| ) for \x — 1| not to large, they 
display a large deviation as \x — 1 increases. This is 
due by the terms neglected in (|43]). We remark, how- 
ever, that these become less and less important as k m 
decreases, as a consequence we expect that the range of 
values of \x — 1 where there is a good agreement with 
( ff(| ) should increase with time. This is indeed observed 
by using data for increasing time in the best fit. Roughly 
the range increases as l/k m (t). 

Finally, we have explored, different types of conserva- 
tion laws represented by Tijk^ with < [i < 2. In all 
these cases the dynamics selects for intermediate times a 
peak at finite values of th wavevector k. 



VI. CONCLUSIONS 

To summarize in this paper we have studied the evolu- 
tion of an TV-component version of the Phase field model 
and shown that the coupling between the massless trans- 
verse modes and the diffusive field produce an instability 
at finite wavelengths. 

Our model, where the low-energy Nambu-Goldstone 
modes couple to the diffusive modes provides an interest- 
ing new scenario and we believe represents a paradigm 
for the Mullins-Sekerka type of phenomena where the 
soft modes are represented by the capillary wave spec- 
trum associated with the solid-melt interface and the 
diffusive mode is the heat transport. These two fields 
concur to destabilize the solid-melt boundary in analogy 
with our findings. On physical grounds, one expects this 
kind of instability to occur during phase separation, be- 
cause small droplets can dissipate heat more efficiently 
and reach rapidly thermal equilibrium due to the Gibbs- 
Thomson effect, whereas larger droplets try to dissipate 
energy faster by creating bulges thus increasing the cur- 
vature. As the system cools down upon reaching equilib- 
rium the typical length of the bulges, A/, increases and 
diverges together with the average domain-size L(t). 

We have demonstrated that the presence of U induces 
non trivial effects on the field cf> because it acts on a time 
scale longer than the noise field, characterized by a short 
correlation time. 

Finally, the O(N) model analysed presents unusual fea- 
tures since it displays scaling behavior in the early regime 
(t < Tf) and multiscaling |l5|] in the late regime and 
constitutes an example of multiscaling without COP, a 



phenomenon which at the best of our knowledge was not 
observed before. 

In summary the present model reveals an unexpectedly 
rich dynamical behavior which reminds in many respects 
the solidification kinetics. 



VII. APPENDIX: EQUILIBRIUM PROPERTIES 

In the present appendix we shall outline the calcula- 
tions of the equilibrium properties of the model. The 
partition function associated with the Hamiltonian 

/oo 
CA^e-^M+W (48) 
-oo 

Where we have included an external field h(x) coupled 
linearly to <fi(x) and (3 — (fc^T/) -1 . The field U(x) can 
be traced out and one finds a part from uninteresting 
constants: 



Z[{h(x)}] 



n„ =1 # Q e 



(49) 



where H c s is defined by Eq. (|25j). In order to separate 
the macroscopic component P of the field we employ the 
following identity 



/oo 
dP 2 5(NP 2 -^2(f>l) 
-OO 



(50) 



and rewrite Z as 

z = n f°° dP 2 f°° & r°° n£_!d0 a 

J— oo J— oo 2ir J— oo a — 1 rot, 

x e -fW ca [<t, a \+l3h<l>+i\(NP 2 -Y^ a <t>l) ( 51 \ 



z=n HdP 2 r ^e-^^+iP*-^ 2 ) 



oo J —oo 



2tt 



/oo 
^ =1 d(p a 
-oo 

xe _/3/2 S«[- 1 d d xV<f> a v<f, a +2i\/p<f> 2 a }+t3J2 a h °'<t'°> (52) 

In the case of a uniform external field directed along the 
component 1 (hi — h), after eliminating the <j) a fields, Z 
reads: 



'"° dX 

2-7T 

oo J — oo 



z = Ne $W**/fi) I dF 2 I ™ e Nn (53) 



Icff P 2 J_ 2. P 4 _ i A P 2 
2 4 C (3^ 



pjd d x 



f) A h A 
4A 



(54) 



In order to evaluate Z we apply the saddle point es- 
timate in the limit N — > oo imposing the conditions 
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(dZ/dX) = and (dZ/dP 2 ) = which lead to the con- 
ditions: 



P 2 



— =r cS + gP 2 



d d k 



(55) 



2 h 2 



P 2 h 



(2ir) d k 2 + 2i\/[3 4A 2 
Eliminating A with the help of eqs. (|5q)-(|5q) we find 



(56) 



P 2 



h 2 



1 



d d k 



1 



(r cff + gP 2 ) 2 (3 J {2it) d k 2 + r cS + gP 2 
m 2 + 5oo 



(57) 



The last term equals m 2 , the square of the average mag- 
netization per unit volume m = ydlnZ^/dh. By using 
eq. (pq) we find explicitly: 



2i\ 



off + gP 2 



(58) 



The existence of a spontaneous magnetic phase implies 
that in zero magnetic external field m ^ 0, i.e. the fol- 
lowing condition must be fulfilled: 
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lim [r eS + gP 2 } = 

h^0 



(59) 



The equation of state reads 

T , d d k 

9 T f 



(2v) d k 2 + r cff + gS 



+r c s + gm 



gm* 



m 



(60) 



FIG. 1. Dispersion relation ui+(k) in the early regime and 
in the late regime. Both quantities are plotted in arbitrary 
units. Notice the different long- wavelength behavior in the 
two cases. The first is typical of NCOP dynamics, whereas 
the second characterizes COP dynamics. 



where Soo is given by: 

Q T f ^ 

S °° - 9Tf J J2^~ d 



1 



k 2 + r cff + gS 



gm A 



(61) 



In order to determine the critical temperature T c we re- 
quire to 2 = and r c ff + gSoc = 0. 



T c = (p 2 /w-r)(d-2)/(gk d - 2 K d ) 



(62) 



where 1/K d = 2n d / 2 T(d/2) where T(x) is the Gamma 
function. 
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FIG. 2. Dependence of the wavevector instability kf on 
the coupling parameter fi in log- log scale. Both quantities are 
plotted in arbitrary units. The squares are obtained from the 
numerical solution of eqs. ( |l8| ) - ( po"[ ) with = 1, Tu = 5, 
r — —0.5, g = 1, w = 0.05 and d — 3. The broken line has 
slope 1. The value of kf is taken when the peak at kf becomes 
larger than the one at k = 0. Other definitions, e.g., fastest 
growing mode, lead to a similar dependence as follows from 
the rescaling (|34j) . 



FIG. 3. Dependence of the crossover time r/ on the cou- 
pling parameter /i in log-log scale. Both quantities are plotted 
in arbitrary units. The broken line has slope —2. The squares 
are obtained from the numerical solution of eqs. ( |l8[ ) - ( pc| ) 
with = 1, Tu = 5, r = —0.5, g = 1, w = 0.05 and d = 3. 
The crossover time 77 is defined as the time when the peak 
at k 7^ becomes dominant. Other definitions, e.g., fastest 
growing mode lead to similar dependence as follows from the 
rescaling ([m]). 
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FIG. 4. Typical evolution of the structure function 
C^(k,t) for different times as shown in the figure. Both 
quantities are plotted in arbitrary units. Notice that the two 
largest times correspond to t > 17 and display the character- 
istic CO P p eak. The data are from the numerical solution 
of eqs. @ - @ with = 1, IV = 5, r = -0.5, 3 = 1, 
w — 0.05 and d — 3. 



FIG. 7. Height of the peak of the structure function 
C^(k m ,t) as a function of time. Both quantities are plotted 
in arbitrary units. The crossover from the NCOP behavior to 
the COP behavior is evident. The dashed line represents the 
t d / 2 behavior whereas the dashed-dot line the t d ^ A behavior. 
The data are from the numerical solution of eqs. jl^ ) - j2c| ) 
with parameters as in Fig. H. 



FIG. 5. Typical evolution of the structure function 
Cd,u(k,t) The data are from the numerical solution of eqs. 
( [L8| ) - (|2y) with parameters as in Fig. ^. Both quantities are 
plotted in arbitrary units. 



FIG. 8. Shape function F{x) of the C^^{k,t) structure 
function in the late stage evolution. The absence of scaling is 
evident in figure. The data are from the numerical solution 
of eqs. ( |l8| ) - ( pp| ) with parameters as in Fig. ^j. 



FIG. 6. Time evolution of the structure function Cuu (k, t) 
for different times. Notice the conserved dynamics at all 
times. The data are from the numerical solution of eqs. (^) - 
( pp| ) with parameters as in Fig. Both quantities are plotted 
in arbitrary units. 



FIG. 9. Multiscaling exponent ip(x) defined in the text. 
The data are from the numerical solution of eqs. ( |l8| ) - ( |2c| ) 
with F^, = 1, Tu = 5, r = —0.5, g = 1, w = 0.05 and d = 3. 
The crosses are obtained from while the triangles from 
Cuu- The full line is the theoretical prediction (|46j) . 
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